Critical Dynamics of Superconductors in the Charged Regime 
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We investigate the finite temperature critical dynamics of three-dimensional superconductors in 
the charged regime, described by a transverse gauge field coupling to the superconducting order pa- 
rameter. Assuming relaxational dynamics for both the order parameter and the gauge fields, within 
a dynamical renormalization group scheme, we find a new dynamic universality class characterized 
by a finite fixed point ratio between the transport co-efficients associated with the order parameter 
and gauge fields, respectively. We find signatures of this universality class in various measurable 
physical quantities, and in the existence of a universal amplitude ratio formed by a combination of 
physical quantities. 



Close to the critical temperature T c of the normal- 
superconductor transition, in a regime determined by 
the Ginzburg criterion 0, 0], order parameter fluctua- 
tions dictate critical properties. For decades, the effect 
of the charge of the superconducting order parameter in 
this regime in three dimensions has formed the topic of 
keen study. While for strongly type-I materials, the cou- 
pling of the order parameter to transverse gauge field 
fluctuations is expected to render the transition first or- 
der 0, it is well-established that strongly type-II mate- 
rials should exhibit a continuous phase transition, and 
that sufficiently close to T c , the charge of the order pa- 
rameter field is relevant Q. While the exact location of 
the boundary between these two types of behavior is still 
the subject of investigation ||, the static critical proper- 
ties of the charged- XY universality class are reasonably 
well-understood. With the discovery of superconducting 
compounds with large critical temperatures and short co- 
herence lengths the critical regime of this transition is 
now potentially accessible to experimental investigation. 

In this Letter, we investigate the less-addressed issue 
0, Q of the dynamics of the three-dimensional normal- 
superconductor transition in the charged regime. In the 
well-studied case of superfluid He 4 , the coupling of the 
order parameter to a conserved energy density field has 
non-trivial effects on critical dynamics Analogously, 
we propose that the coupling of the superconducting or- 
der parameter to relaxational transverse gauge field fluc- 
tuations leads to qualitatively new dynamics character- 
ized by a universal ratio C between the zero wave-vector 
part of the characteristic frequencies for the dynamics of 
the order parameter t/j a and that of the gauge field A at 
the critical point: 



const. 



(1) 



Thus we propose the possibility that the strong coupling 
of these two fields causes them to relax in the same fash- 
ion at the critical point, with a single new dynamic ex- 
ponent z. In what follows, employing dynamic renormal- 



ization group (RG) techniques, we show these features to 
hold within the context of a particular model. We discuss 
the universal properties obtained from the model and the 
behavior of measurable physical quantities. 

The Model- As a starting point for modeling dynam- 
ics, we employ the standard Ginzburg-Landau free energy 
used to define the finite temperature static critical prop- 
erties of a three-dimensional superconductor coupled to 
a transverse electromagnetic field 0, generalized to N 
complex species of matter field, given by 



F = I d d x 



V x A| 2 



(2) 



where ipa (with a G 1,2, ...A) is the generalized order 
parameter and A is a fluctuating massless gauge field. 
The effective charge is given by Jg = e*\/4Tr/hc and 
deviations from criticality are measured by r = 2m* (T — 
T c )/h 2 T c [13]. For N = 1 and d = 3, this model describes 
the low-energy excitations of a bulk charged superfluid in 
the regime where charge-density fluctuations are gapped 
at a high-energy (the plasmon gap), but remain coupled 
to transverse electromagnetic fluctuations. 

As originally obtained in Ref. || , a one-loop static RG 
analysis of Eq.JSJ) in d = 4 — e dimensions shows that 
no non-trivial charged fixed point exists for 2N < n c — 
365.9. However, for a range of parameters, more sophisti- 
cated methods indeed find a continuous phase transition 
for N = 1 [3. IToL FTl| - describing, presumably, the second 
order transition found in many materials. Many of the 
salient features of this "charged- XY" universality class 
are captured by Eq. (J2Jl with N > N c ~ 183, for instance 
the fact that the anomalous dimension of the order pa- 
rameter field is negative. Given that the AT-component 
model does provide insight on the statics of the charged 
superconducting transition, here we make the reasonable 
assumption that when the free-energy defined by Eq.J2J 
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is augmented by appropriate equations of motion, the 
model captures the basic features of the critical dynam- 
ics of the transition as well. 

For completeness, we note that the free energy density 
is invariant under the following transformations: 



1pa — > SabOapTpp, 

A(r) -> A(r)+VA(r) 



(3) 
(4) 



where a, (3 e {1, ...,N}, a,b E {1,2}, O is an N X N 
orthogonal matrix, 9 E [0, 2tt], and a v ab is the usual 2x2 
Pauli matrix, rp a and tp^ are the real and imaginary 
parts, respectively, of the a-th complex field. The term 
in Eq.® with coupling constant u is the only quartic 
matter-matter interaction allowed by these symmetries. 

The simplest equations of motion augmenting the stat- 
ics described by Eq.® are relaxational in both ip a and 
A. In the presence of external fields h^ a and hj^ coupling 
to ip a and A respectively, they take the form 

d t tpa = -I> (J^ - +?7 Q , (5) 

d t A t = -T A (J^- - h A ^j + Q. (6) 

Here, and Ta are transport co-efficients associated 
with ip a and A respectively. The fields r] a , ??* and £ are 
white noise correlated, and ensure that the fluctuation- 
dissipation theorem is satisfied. Thus, they obey the con- 
straints 

(?j a (r,t)^(r',t')) = 2TT^6 aP 6 3 (r-r')S(t-t'), (7) 
fa a (r,t)) = (r, a {v,t)r] (r',t')) = 0, (8) 
(CiM)g(r',t')) = 2TT A 8 t3 8 3 (r-r')5(t-t'), (9) 
(Ci(r,*)> = 0. (10) 

In fact, if we assume that the normal state of the system 
is a metal, Eq. ® can be derived from the low- frequency 
form of Maxwell's equation 



V x B = 4ttJ/c + dtE/c, 



(11) 



in the gauge Aq = 0, where the electric and magnetic 
fields are given by E = — 9 t A/c and B = V x A, re- 
spectively. The net current J = j s + j„ has a super- 
fluid component j s = ~ iy/g-^^a] 2 and a nor- 
mal component j n . The average normal current is given 
by <t„E, where a n is the normal conductivity. Thermal 
fluctuations of the normal fluid give rise to the noise in 
Eq.©. With these assumptions, one can retrieve Eq.® 
from Ea. (|llf> in the limit u> — > 0, and identify the inverse 
transport co-efficient, T^ 1 , with the bare conductivity, 

On- 

A complete model for the dynamics requires identifi- 
cation of all conserved quantities and Poisson-bracket re- 
lations applicable to the normal-superconductor system 



in the charged regime [8J. Even in the uncharged regime, 
one might expect non-dissipative coupling of the order 
parameter to a combination of energy and mass densit 
to exhibit model E dynamics as in superfluid helium 
However, hydrodynamic analyses show that in the pres- 
ence of impurities (which is implicit in the assumption of 
finite conductivity well within the normal state) this cou- 
pling does not survive [12|, as indicated by the absence 
of second sound modes in actual superconducting sys- 
tems. However, in principle, a conserved energy density 
mode could couple to the order parameter mode via non- 
linear interactions, leading to model C dynamics in the 
uncharged superconductor [l^ . Likewise, in the charged 
regime, we do not expect any non-dissipative coupling of 
energy-mass density to the order parameter, due to the 
presence of impurities in real samples of interest. How- 
ever, the possibility of non-linear interactions with con- 
served quantities cannot be completely eliminated. The 
model we employ consisting of the order parameter and 
gauge field modes alone, each with relaxational dynam- 
ics, is the simplest, but hitherto unexplored, possibility. 

RG analysis- The equations of motion Eqns.(0EJ 
allow a dynamical RG analysis which we now detail (also 
see, for example, Ref.||). The effective charge ^fg and 
the coupling constant u of Eq.® are treated perturba- 
tively, as is the deviation from four dimensions, e = 4 — d 
in order to avoid infrared divergences 0|. Because the 
characteristic electron speeds are small compared to the 
speed of light, the low-energy theory of the system need 
not be relativistically invariant. This leads to the resid- 
ual gauge symmetry contained in Eq. . (This has been 
pointed out previously in, for instance, Ref. 14].) To 
avoid divergences in functional integrals resulting from 
the fact that multiple choices of gauge lead to the same 
magnetic field configuration, we perform the Fadeev- 
Popov procedure and add a term 1/(2A)(V • A) 2 to the 
free energy (Eq.®) 15]. Fixing the value of A corre- 
sponds to making a choice of gauge, and for the remain- 
der of the paper we work with A = 1 (analogous to the 
Feynman gauge in quantum electrodynamics). Regard- 
ing the dynamics, since our focus is on the manner in 
which the relaxational rates of the fields ip and A af- 
fect one another, we rescale the theory and write the 
equations of motion in terms of the ratio of transport 
co-efficients T = T^/Ta- 

We perform the standard RG procedure, integrating 
out modes in a momentum shell A/6 < |k| < A, where A 
is a high momentum cut-off, and all frequencies, followed 
by a rescaling of space and time: r — > br and t — > b z t. To 
one-loop, and 0(e), the j3- functions for u and g are (as 
in Ref. 0]): 



du 
dlnb 

dg 
dlnfe 



eu - 2 (AT + A)u 2 - 3g 2 + Qgu, (12) 
N 



(13) 
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yielding Gaussian (u* = g* = 0), XY (u* ^ 0, g* = 0), 
and charged- XY (u* ^ 0, g* ^ 0) fixed points for N > 
183. At the charged-XY fixed point of interest, g* — 
3e/N to leading order in e. We find that the anomalous 
dimensions of the order parameter and gauge field at this 
order (and in our choice of gauge) are: 



i.e. 



-2g ; VA = Ng/Z, 



(14) 



and note that t\a = e = 4 — d at the critical point, as 
required by gauge invariance ^lj • 

We fix the exponent z by requiring that the time 
derivative term in the equation of motion for A return to 
its "bare" form. The same feat in the equation of motion 
for %/} is accomplished by allowing the ratio T to flow. We 
find that 



1 = tf-^-zfl+g^bY 



r' = b z - 2+,, *-r i 
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giving a one- loop (3- function for T: 
dT / N r 



r + i 



In b 



d\nb 



N 



2v r + i 



(15) 
(16) 

(17) 



At the 0(e) charged fixed point, Eq. l|17|) has a stable 
fixed point solution: 



(JV- 12) + V(A - 12)^ + 2AN(N + 9) 

r W+v) ' ( } 

which in general obeys < T* < 1.5 for all N > 0. 
We note that T = (T^ — > with fixed, for in- 
stance) is not a fixed point and that = (r^ — > 
with r^, fixed, for instance) is an unstable fixed point, at 
this order. Thus, we are led to conclude that, at least 
within a one- loop RG analysis, the critical dynamics of 
the charged superconductor is governed by a non-trivial 
fixed point wherein T has a finite ratio, reflecting the fact 
that the dynamics of the order parameter and those of 
the gauge field are strongly coupled. At the charged-XY 
fixed point, we find that z = 2 + e[3/(2T*) - 1] to lead- 
ing order in e, implying that z > 2. While an accurate 
value of the dynamic exponent would require employing 
more sophisticated treatments, our results certainly sug- 
gest that near criticality, the system relaxes slower than 
expected for diffusive dynamics. 

Physical Consequences- The most striking new fea- 
ture of this fixed point (in contrast with the uncharged 
dynamics) is the existence of the universal fixed ratio, T*. 
Physically, its existence requires that the order parame- 
ter and gauge fields relax in the same fashion. In fact, 
using scaling arguments one can show that at the 
critical point, the ratio of the characteristic frequencies 
at k = of the two fields is exactly the universal ampli- 
tude ratio, identifying the constant C of Eq.QJ with T*, 



C^hm^TT 



(19) 



Here, the characteristic frequency of a field Q is defined 
by 



w Q 1 (fc) = iXQ(k,u = 0) 



(20) 



where XQ is th e dynamic linear response function, (Q) = 
XQ^Q- I n our case, close to criticality, the characteristic 
frequency has the scaling behavior uiq = ^ z flQ(k£ i ) for 
both fields "0 and A, where £ is the divergent correlation 
length associated with the order parameter, and Jig is 
a universal function. In principle, the characteristic fre- 
quency for each field can be obtained by measuring the 
static susceptibility and dynamic linear response function 
associated with the field. 

The dynamic linear response functions themselves 
carry valuable information on critical dynamics. Close 
to T c , they have the scaling form 



X Q (k,uj)^e- riQ fQ(^ z ,kO, 



(21) 



where, associated with each field ip a and A, t]q and /q 
are the anomalous dimension and a scaling function, re- 
spectively 

Each of these functions is manifest in measurable quan- 
tities. The order parameter response function, Xip(k, w), 
is the pair susceptibility appearing in Josephson tunnel- 
ing experiments |l6| . In principle, the scaling behavior 
of Xip m Ea. (|2"T|l could be used to extract the dynamic 
exponent, z. Accessing xa(^ = 0, us) should be relatively 
straightforward, since it is related to the resistivity in 
linear response. The external field that couples to A 
is an applied current j ex t, so that, (A) = XAjext, with 
a(k = 0, u>) = — ix^ 1 (k = 0,w)/w. Close to criticality 
and above T c , we expect this to provide the dominant 
contribution to the net conductivity and, using Ea. H21|) . 
to diverge as a{k = 0,lu) = £, z ~ 2+r,A G{u£, z ) - £ 2+2 - d , 
where Q is a scaling function. This is a consequence of 
the Josephson scaling relation 0, ^| which holds in 
both charged and uncharged regimes, as a result of gauge 
invariance 11]. Notably, in the charged regime, the ex- 
ponent z has a different value than in the uncharged case. 

We see that various features of the dynamic universal- 
ity class of the charged superconductor appear in mea- 
surable quantities. The definitive signature of this uni- 
versality class would be in the extraction of the universal 
amplitude ratio, Eq. l|19|) . 

Concluding Remarks- The charged 

superconductor-normal transition is a continuing 
source of rich physics. Research on the static transi- 
tion of the charged Ginzburg-Landau model indicates 
differing behavior for type-I and type-II superconduc- 
tors as well as the relevance (sufficiently close to the 
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continuous transition) of the charge of magnetic field 
fluctuations. Numerical and analytical work in Ref. 6] 
on the dynamics of this transition also seem to reveal 
distinctly different dynamic behavior between strongly 
and weakly screened superconductors. In Ref.0, Monte 
Carlo studies of superconductor dynamics in the vortex 
representation find that for strong magnetic screening, 
z rs 2.7, in qualitative agreement with our result that 
the order parameter dynamics in the charged regime is 
sub-diffusive. Here, as our key point, we suggest that for 
materials with continuous transitions from normal metal 
to superconductor, the dynamics in the charged regime 
will be governed by a new universality class. Further 
analyses of all these issues are well in order. 

In experiments, the charged regime of the 
superconductor-metal transition is not easily acces- 
sible. The Ginzburg criterion indicates that materials 
with high critical temperature, large anisotropy and 
extreme type-II behavior should manifest large regimes 
of fluctuations. However, within this fluctuation regime, 
the region close to T c where the system crosses over to 
the regime of charged fluctuations is often too narrow 
to access Q. For instance, high T c materials such as 
YBa 2 Cu307_(j, while possessing large regimes of critical 
fluctuations, are too strongly type-II to observe charged 
critical fluctuations. However, weakly type-II materials 
with high T c 's or granular texture, and moderate 
anisotropics could open up a window onto this new 
regime. We are hopeful that an investigation of such 
materials will yield an understanding of the effect of 
charge on the critical dynamics of this transition. 
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